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Abstract. We explain why the macroscopic fluctuations of deterministic local 
collision dynamics should be characterized by a non strictly convex functional. 

The study of thermal conduction properties of extended lattice Hamiltonian sys- 
tems has received recently a lot of attention by theoretical and mathematical physi- 
cists alike. The first issue is simply the derivation of the Fourier law, or equivalently 
the heat equation. On the other hand, an impressive amount of rigorous results 
concerning stochastic interacting particles system has been obtained over the last 
decades. The diffusion equation has been derived in the hydrodynamic limit and the 
behaviour of the fluctuations of the macroscopic observables of the system under a 
space-time diffusive scaling has been derived rigorously . In particular, the shape of 
the large deviations functional of the current and density of particles profile has been 
determined [H El El IU El El [3, [Hj. This object may be seen as the out-of-equilibrium 
analogue of the equilibrium thermodynamic potentials. 

The first example of local collision dynamics has been introduced in [8j. In that 
model, particles are locked in cells having a particular shape ensuring strong chaotic 
properties for the dynamics of each particle within its own cell. The cells are ar- 
ranged so as to tile the plane and have a small opening so that particles in neighbour- 
ing cells may interact through elastic collisions. This was model taken up in [5] as 
a model for heat conduction when its boundaries are connected to heat baths. The 
authors argue that in a weakly interacting regime and after a suitable time rescaling, 
the energy exchange between neighbouring cells is described by a stochastic dynam- 
ics analogous to the ones encountered in stochastic interacting particles. This raised 
the hope of using the rigorous results obtained for those systems to derive Fourier 
law for deterministic dynamics. Unfortunately, the stochastic dynamics describing 
the exchange of energy does not satisfy the gradient condition and is therefore diffi- 
cult to analyze. In this paper, we want to explain briefly how another approach to 
local collision dynamics allows to conjecture important features regarding the fluc- 
tuations (large deviations) of its macroscopic observables. This was introduced and 
developed in [101 EED, [El US]. Consider N particles of unit mass with positions and 

momenta (q, p) = {(qj, Pi)} 1<i<N , with q,,Pi G M. d . The positions are measured 




with respect to iV fixed centers located on a ID lattice. The Hamiltonian H takes 
the form : 



where the interaction potential U is equal to zero inside a region VL V C R d with 
smooth boundary A of dimension d—1, and equal to infinity outside. Likewise, the 
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pinning potential V is assumed to be zero inside a bounded region Qy and infinity 
outside, implying that the motion of a single particle remains confined for all times. 
The regions Qjj and Qy being specified, the dynamics is equivalent to a billiard in 
high dimension. A typical example of the dynamics we wish to consider is given 
by the figure below. The circles move freely within their square cells and collide 
with each other when they both get sufficiently close to the hole located in the wall 
separating two adjacent cells. 




Physically, those models describe aerogels, i.e. gels whose liquid components have 
been removed and replaced by atoms of gases. One of the simplest example of this 
type of dynamics is the complete exchange model introduced in [15]. Although it 
is expected to display anomalous thermal conduction properties, it has nice aspects 
that underline the general important features of local collision dynamics. In the 



definition of the Hamiltonian (0.1), take d — 1 and the potentials V and U given by 



v(r) _ f +oo if \x\ > b _ f +oo if \x\ > a 

V{X) ~ \ if \x\ < b U{X) ~ \ Oif \x\ < a 

Each particle on the lattice moves freely on a one-dimensional cell of size 2b, 
changing directions at the boundaries. The interaction between a pair of particles 
acts when the difference between the positions of the two particles reaches the value 
a, at which point they exchange their velocities. 

The motion of a given pair of particles at sites i, % + 1 is described as the motion of 
a point particle on a two-dimensional billiard table described by 

f2 = {(xi, Xz) G R 2 , \x\\ < b, |x 2 | < b, \xi — x%\ < a} . 



26 




Figure 2. Billiard giving the motion of two particles in the complete 
exchange model 



It is straightforward to see that the energy of the particle n at time t is given by 
E(n,t) -E(n,0) = J(n-1, [0,t]) - J{n,([0,t)) 
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where J(n — 1, [0, £]) is the time-integrated current: 

, c n (t) 

J(n, [0,t]) = --^[^ +1 (r n fc )-^(r„ fc )] 
fc=i 

(r^)fc is the sequence of collision times between particles n and n + 1 and C n {t) is 
the number of collisions up to time t. For more general local collision dynamics, 
the current keeps this "gradient" -like expression with a prefactor in front of the 
difference of kinetic energy, see (10] for a detailed expression. This a fundamental 
difference with the expression of the current that is found in chains of oscillators 
interacting through "soft" interactions (i.e interacting with a smooth potential U). 
In those systems, the time integrated current is of the form 

1 /•* 

J(n, [0,t]) = -- J (p n (s) + p n+ i(s))U (q n (s) - q n+1 (s))ds, 

when the particles move in a Id potential wall. It is easy to see that in that case, 
the expectation of the instantaneous current (the integrand in the above expression) 
with respect to a local equilibrium measure is zero. This is not the case with the 



expression in (0.2). 



We explain now how to construct a stochastic dynamics that is simple to ana- 
lyze and yet hopefully retains the essential features of the determinisitic models. 
The basic idea behind the model is the following : when the system is large, each 
particle in a given cell sees its neighbouring particles as being embedded in a heat 
bath whose temperature changes slowly on a microscopic time scale, meaning that 
many collisions occur before the temperature of the bath is significantly affected. 
Conceptually, one should note that it is a very different way (and a more generic 
one) of introducing randomness in the picture than the one advocated in [9]. Here 
the randomness in the motion of one given particle appears as the result of the 
action of a large number of particles, not as an effect of the specific shape of the 
cell in which each particle moves. Roughly speaking, we replace the walls through 
which the collisions between neighbours occur by stochastic thermal walls at fixed 
temperatures and let the particles travel back and forth between the thermal walls. 
In fact, we dissociate two aspects of the dynamics : each particle carries some energy 
by travelling between the two sides of its cells and on the other hand, it acts on its 
neighbours as part of a huge thermal bath when colliding with them. In summary, 
we build a model made of tracers and scatterers, the latter modelling the actions of 
a heat bath at a given temperature. Pictorially, we get this : 
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Figure 3. The motion of a particle between two thermal walls 
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We delete the vertical dimension that is, for our purpose, irrelevant. And we first 
concentrate on the motion of one particle between two thermal walls located at the 
ends of the interval [0, 1]. We denote its position and velocity by (q(s),p(s)). 



Figure 4. The motion of a particle between two thermal walls 



The motion is ballistic inside the interval and when the particle encounters one 
wall, its velocity is reversed and its speed is updated with the law : 

<P±(v) = Pve'^ , ± = T ± . 
One may associate to the motion of the particle the time-integrated current : 



2 ^2 

k=l k=l 



„2 



where ~ V 9± ( t ') = /3±ve /3±V ' 2 and JV^ counts the number of visits to the ±-side 
of the interval. Regarding this current, one can show the following : 

J[0,t] T_-T + 
hm — — - = ; r a.s. 

where T_ and T + are the left and right temperatures. But one can go one step 
further and study the fluctuations of the LHS before the limit is taken. Namely, 
one may compute the large deviations functional of the current r, T) roughly 
defined as, 

P r t is the stochastic dynamics with a fixed temperature difference r = T_ — T + and 

2 



average temperature T — T +T+ 



If r 7^ then one gets the following scaling result, 

j — KT 



( (j-KT) 2 -r ■ 2 
a kT 2 11 ]T > KT 



hme- 2 I(ej,er,T) = g(j,T,T) 



if jr e [0,kt 2 ] 
—£s if jr e [-/tr 2 ,0] 



2 

4kT2 11 J 1 ^ ~ KT ' 



where k = (^) 2 - The reason for the affine part of the functional is the fact that 
the particle will get a small velocity with a not so small probability and this will 



MACROSCOPIC FLUCTUATION THEORY OF LOCAL COLLISIONAL DYNAMICS 5 



. 2.0 
\ 15 




\ 1.0 




\ 0.5 





-3-2-10 1 2 3 



Figure 5. Plot of Q clS db function of j for kt = kT = 1 

make the occurence of small values of N t (and thus of J([0, £])) quite likely. This 
explains the fact that the functional vanishes between and the average value kt 
: it is not exponentially unlikely to observe values of the current smaller than its 
average. The scaling expressed in the relation above turns out to be exactly what we 
need to conjecture the shape of the large deviations functional in the original billiard 
dynamics. Indeed, we want to mimick the slow variation of the local temperatures 
in the deterministic model, i.e it is expected to change by a significant amount 
at a speed ~ 1/N 2 at the microscopic time scale. In order to do so, we let the 
temperatures {T(n, t)}o<n<7v of the scatterers in figure 3 evolve such that they satisfy 
the equation 

T(n, t + N 2 At) - T(n, t) = J(n, [t, t + N 2 At}) - J(n - 1, [t, t + N 2 At}), 

where the currents entering the equation are the ones on the LHS and RHS of the 
scatterer n. Define T N : [0, l]xl+4 M+ and J N : [0, l]xl+4 K. 

T N (x,t) = T([Nx],N 2 t) 

J N {x,t) = N-j^- t -J(n,[N 2 t,N 2 (t + At)]), 

Then, one can show[] : when iV — > oo, At — > 0, {Tn-,Jn) converge in L 2 to the 
unique solution (T, j) of 

f d t f(x,t) = -d x J(x,t) 

\ J(x,t) = -K(f(x,t))d x f(x,t) 

with k(T) = (^-)» and suitable b.c. 

At finite N, and for each (x,t) G [0, 1] x [0, 1], J~]sr(x,i) and T^(x,t) are random 
variables and the object one is really interested in is the large deviations functional 
X appearing in : 

F({T N ~f,J N ~j}on[0,l] x [0,1]) ~exp[-iVX(j,f)]. 



Precise statements as well as rigorous proofs of the results below will be the subject of an 
upcoming publication 
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and one can show that 

i(j,f)= [ dt [ dxg(j(x,t),d x f(x,t),f(x,t)) 

Jo Jo 

if j and T satisfy d s T(x, s) = —d x j(x, s), and I = +00 otherwise. The integrand Q 
is given by, 

f ^ if JT > KT 2 

if jt e [0, KT 2 ] 



if JT e [~KT 2 ,0] 



if j T< _ KT 2 



4kT 2 



And therefore, 

which shows that local billiard dynamics should have macroscopic fluctuations prop- 
erties very different from the ones observed in the usual stochastic interacting par- 
ticles such as the exclusion process. The large deviation functional is not a strictly 
convex function of the current. This aspect should translate into some interesting 
phase transitions-like effects. 

The stochastic dynamics that is introduced is a very crude approximation of the 
deterministic dynamics. However, one should note that it is not the fact that we 
replaced the deterministic collisions by stochastic ones that is responsible for the 
special behaviour of Q. This is caused by what is left of the deterministic dynamics, 
namely the ballistic motion of the particle and the fact that one assumes the particles 
to be in local equilibrium. Therefore, one should expect this feature to be robust 
and independent of the approximations that were made for the collisional part of 
the dynamics. 
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